Core excitation in the elastic scattering and breakup of n Be on protons 
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The elastic scattering and breakup of n Be from a proton target at intermediate energies is studied. 
We explore the role of core excitation in the reaction mechanism. Comparison with the data suggests 
that there is still missing physics in the description. 
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I. INTRODUCTION 

Nuclear reactions offer the most diverse methods to 
study nuclei at the limit of stability. Understanding re- 
action mechanisms in nuclear processes involving nuclei 
near the driplines is of great importance, particularly at 
this time, when there is such a high demand for accuracy 
on the structure information to be extracted from the 
data. Reaction and structure models are undoubtedly 
entangled, therefore improving reaction models often im- 
plies incorporating more detailed structure models in the 
description [1J. 

It is generally accepted that, in reactions with loosely 
bound nuclei, the coupling to the continuum needs to be 
considered. Continuum effects are very much enhanced 
in breakup but can also have imprints on other reac- 
tion channels, for example elastic and inelastic scatter- 
ing. One framework that explicitly includes continuum 
effects is the Continuum Discretized Coupled Channel 
(CDCC) method Q. A large amount of work has been 
devoted to the analysis of experiments within this frame- 
work @, 0, [E IS 0] an d m general results are very good. 

The extended Continuum Discretized Coupled Chan- 
nel (XCDCC) method @, Q was recently developed. It 
brings together a coupled channel description of the pro- 
jectile with a coupled channel model of the reaction, en- 
abling the description of interference between the multi- 
channel components of the projectile as well as dynam- 
ical excitation of the core within the projectile, during 
the reaction [9|. The model has been applied to the 
breakup of n Bc^ 10 Be-|-n and 17 C^ 16 C+n on 9 Be at 
re 60 MeV/nucleon gQ. The projectiles are described 
within a two-body core + n multi-channel model, where 
the core can be in the ground state but also in an ex- 
cited state. This model produces breakup cross sections 
to specific final states of the core, given a coupled chan- 
nel Hamiltonian for the projectile. Results presented in 
Ref. [8|, l9j show that core excitation effects in the total 
cross section to the ground state of the core are small, 
but become very large when considering the total popula- 
tion of the core's excited state. Other differences can be 
seen in angular and energy distributions but at present 
no such data is available. The effect of core excitation 



needs to be studied in other regimes, for very light and 
very heavy targets, as well as a variety of energies. 

In this work we concentrate on the proton target. In 
this case the process is nuclear driven, and recoil effects 
are very important. Several reactions of 11 Be on protons 
have been measured in a number of facilities, namely 
elastic scattering at 50 MeV/nucleon [l(|, quasi-elastic 
and breakup at 64 MeV/nucleon [ll[ as well as transfer 
at 35 MeV/nucleon [Hj]. 

n Be proton elastic scattering at 49.3 MeV/nucleon 
was performed in GANIL at the same time as the 
elastic scattering of the core 10 Be at 59.4 MeV/nucleon. 
Even though the outcoming 11 Be measurements corre- 
spond to quasi-elastic, these are essentially elastic as 
the contribution from the first excited state is negligible. 
Standard optical potentials (either density folding as in 
JLM [l3[ or global optical potentials coming from elas- 
tic fits as in CH89 [1J|) could reproduce the 10 Be elastic 
reasonably well, requiring small renormalizations of the 
real and imaginary parts of the interaction (Ay = 0.9 and 
Xw — 1-1 for CH89) [10( ■ Larger renormalizations were 
required in order to reproduce the distribution of 11 Be 
(Ay = 0.7 or X w = 1.3 for CH89) [To|. 

It is clear from Ref. [T(| that the global optical poten- 
tial overestimates the elastic cross section for 11 Be. In 
Ref. [IH the elastic scattering of 11 Be on 12 C was success- 
fully described using a 10 Be+n two-body model, incor- 
porating breakup effects. As the 10 Be-target interaction 
was fixed by the 10 Be elastic scattering data, the large 
modification in the 11 Be+ 12 C elastic data was described, 
without rcnormalization, purely through breakup effects. 
Due to the loosely bound nature of the last neutron in 
11 Be this loss of flux from the elastic channel can be at- 
tributed to breakup into 10 Be+n. It is thus possible that 
the large renormalizations for 11 Be scattering on protons 
[l(| are also due to breakup effects. 

In Ref. [ll|, elastic data is only described after large 
renormalizations of both the real and imaginary part of 
the 10 Be+p interaction (Ay = 0.75 and Xw = 1-8), much 
larger than those used in Ref. [l(| • These same renormal- 
izations can no longer describe the 10 Be+p elastic data 
from Ref. [l(|, and are inconsistent with few-body reac- 
tion theory. We will re-examine the elastic scattering of 



2 



10 / n Be+p to see if one can consistently describe both 
sets of data using the same interaction for 10 Be+p, by 
including continuum and core excitation. 

In addition to elastic and inelastic measurements, 
breakup data from NSCL exist at 63.7 MeV/nucleon 
This breakup data is integrated into two wide energy 
bins due to statistics. The lower energy bin covers the 
1.78 MeV resonance and a reasonable angular distribu- 
tion is obtained, which underpredicts the cross section 
[ill ]. The higher energy bin covers resonances that are 
thought to be built on excited core states. The cal- 
culations presented in Ref. [llj failed to reproduce the 
shape of this higher energy bin, and the authors suggest 
that the source of the disagreement may be due to an 
active core during the reaction. Now that it is possi- 
ble to include core excitation in the reaction mechanism 
Q we will re-examine the breakup data using a consis- 
tent 10 Be+p interaction, and including systematically the 
coupling to the 2 + state in 10 Be. 

Transfer reactions have also been performed with the 
n Be beam, at 35.3 MeV/u in GANIL [3 with the aim 
of extracting spectroscopic factors for the ground state. 
While the reaction mechanisms proved to be more com- 
plicated than the 1-step DWBA theory, results for (p,d) 
show evidence for a significant core excited component. 
The inverse reaction, 10 Be(d,p) 11 Be, has also been stud- 
ied in GANIL [lj], the main interest being the resonance 
structure of 11 Be. This illustrates how transfer is being 
used beyond the standard application of spectroscopy of 
bound states, underlining the need to better understand 
the transfer mechanism and its coupling to the contin- 
uum. 

All these different data offer a good testingground for 
theory. A comprehensive theoretical study [17| focusing 
on 10 Be(d,p) show inconsistencies of the extracted spec- 
troscopic factors for data at different energies. Optical 
potential uncertainties and core excitation effects could 
be at the heart of the problem. 

In this work we perform calculations including elastic, 
inelastic, and breakup channels of 11 Be on protons at in- 
termediate energies. We explore explicitly the effect of 
the inclusion of core excitation in the reaction mecha- 
nism. Comparison to elastic and breakup data will be 
presented here. The analysis of the inelastic channel is 
presented in [l8[ and we leave a detailed study of the 
transfer channel for a future publication. In section |TT] 
we provide the details of the calculations. In section [TTTl 
we present the results: first for the clastic channel pil A|) . 
then for the breakup ((III Bp . Finally, in section IIV1 we 
draw our conclusions and provide an outlook into the 
future. 



energy 


V Rv av W Rw aw 


40 
60 


60.84 1.000 0.7 23.16 0.600 0.6 
31.64 1.145 0.69 8.78 1.134 0.69 



II. 



DETAILS OF THE CALCULATIONS 



TABLE I: 10 Be-proton Woods-S axon potential parameters. 
All energies are in MeV and lengths in fin. 



The calculations for breakup of loosely bound systems 
on protons have a rather different convergence require- 
ment as compared to the breakup on heavier systems. 



The model space needs to span large excitation energies, 
while the radial dependence can be reduced significantly. 

For the CDCC calculations at 40 MeV/nucleon, the 
continuum was discretized upto 35 MeV, with 10 bins 
upto 10 MeV for s-, p-, and d- waves, and 8 bins from 
10-35 MeV. We include 12 bins from 0-35 MeV for all 
other partial waves up to Z max = 4. The same binning 
scheme was used for the XCDCC calculations, except 
that the higher bin density upto 10 MeV was only used 
for channels with outgoing ground state core components. 
Partial waves up to / max = 4 were used for the coupled 
channels projectile states. 

For the 60 MeV/nucleon CDCC calculations, a slightly 
different binning scheme was adopted to match the ex- 
perimental energy bin integrations. From 0-0.5 MeV, 2 
bins were used; over the observed energy bins from 0.5- 
3.0 and 3.0-5.5 MeV, 3 bins were used in each case; and 
from 5.5-30 MeV, 6 bins. For the XCDCC calculations 
where the outgoing channel had excited core states, only 
1 bin was used from 0-0.5 MeV, 1 bin for each observed 
energy range, and 5 bins above. 

The radial integrals for the bins were calculated upto 
40 fm in steps of 0.1 fm. The radial equations in the 
CDCC method were calculated for 30 partial waves with 
the lower radial cutoff for the integrals set to 4 fm inside 
the point Coulomb radius, and matched to the asymp- 
totic Coulomb functions at 150 fm. 

The 11 Be bound state potential parameters are taken 
from Ref. [19], using the Bel2-pure interaction for the 
CDCC calculations and the Bel2-b for the XCDCC cal- 
culations. The 10 Be-proton interaction is fitted to the 
proton elastic data available at the two energies. A good 
fit could be obtained from a renormalized CH89 interac- 
tion [l4j. The parameters are given in Table HI For the 
cases including 10 Be excitation, the OM potentials were 
deformed with the same /?2 deformations as used in the 
n Be bound state. The coupling matrix elements to the 
excited state in 10 Be assume a rotational model with the 
deformation fitted to the experimental B(E2) strength 
[20| . The deformation length is in good agreement with 
that obtained from inelastic scattering ot the 2 + state in 
10 Be [2l| , and the optical potential used here reproduces 
the angular distribution of the inelastic scattering well. 
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FIG. 1: (Color online) 10//11 Be elastic scattering on a proton 
target at ~40 MeV/nucleon, using an optical model fit to the 
10 Be elastic and various 11 Be reaction models for the n Be 
data. The experimental data are from GANIL [22l ]. 
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contributions. 

Calculations were repeated for a higher energy, around 
60 MeV/nucleon, where both 10 Be and 11 Be elastic data 
exist. Once again, when the 10 Be data is fitted with an 
optical model, and the 11 Be elastic is described within the 
CDCC approach, the cross section is over-estimated (see 
Fig. [2]) . Note that the data at this energy does not span 
a large angular range, but it is evident that the pattern 
of over-predicting the 11 Be cross section remains. 

Other reaction calculations have been performed in an 
effort to describe this data [23|, which consisted of a 
transfer to the continuum approach in which the breakup 
continuum was described using the deuteron basis. This 
also failed to describe the data when the 10 Be potential 
was fixed to the elastic data. The same pattern of over- 
predicting the 11 Be elastic was also seen at a lower energy 
of -40 MeV/nucleon [H|. 

As pointed out earlier, in pxl large renormaliza- 
tion factors were needed to reproduce the elastic cross 
section. By including more relevant reaction channels, 
one might account for a part of the renormalization re- 
quired, corresponding to flux that is being removed from 
the elastic channel. This suggests that there are still 
channels coupled to the elastic that have not been con- 
sidered. Preliminary calculations including the deuteron 
transfer channel along with the breakup in the 11 Be basis 
show improvement at small angles, but the disagreement 
still remains at large angles. Due to large non-orthonality 
corrections, CDCC calculations including the deuteron 
transfer coupling turn out to be numerically challenging. 
They will be discussed in a later publication. 



FIG. 2: (Color online) 10 ^ 11 Be elastic scattering on a proton 
target at ~60 MeV/nucleon, using an optical model fit to the 
10 Be elastic and various n Be reaction models for the n Be 
data. The experimental data are from Ref. [hj ( 10 Be) and 
Ref. H3 ( n Be). 



III. RESULTS 

A. Elastic channel 

The elastic scattering is the first test on the reaction 
model. In Fig. [1] we show the 10 Be and 11 Be elastic data 
and theoretical calculations at ~40 MeV/nucleon. The 
optical model for the 10 Be (dashed/black line) is fitted to 
the 10 Be elastic data (open circles). The cluster folding 
model (dotted/red line) folds the 10 Bc+p and n-p interac- 
tions over the 11 Be ground state wave function to produce 
the 11 Be+p potential. Also shown in Fig.[T]is the effect of 
the 11 Be continuum within CDCC (dot-dashed/blue line) 
and core excitation within XCDCC (solid line). Even 
though there is significant improvement over the simple 
optical model when including breakup, results still over- 
estimate the 11 Be elastic cross section at larger angles, 
and no improvement is found by including excited core 



B. Breakup: comparison with data at ~60 
MeV / nucleon 

Breakup data was also obtained at 63.7 MeV/nucleon 
PH . summed into two energy bins. The first covers the 
energy range 0.5-3.0 MeV, which spans the 1.78 MeV res- 
onance, predominantly a d-w&ve neutron coupled to the 
ground state of the core. The second energy bin is over 
the energy range 3.0-5.5 MeV, which spans a resonance 
at 3.89 MeV, thought to be predominantly an s-wave 
neutron coupled to a 10 Be(2+) core [lj. In Ref. p~0], 
CDCC results were presented which underestimated the 
cross section for the lower energy bin, but did not repro- 
duce the higher energy bin. It was suggested that since 
the higher energy bin spanned a resonance with a pos- 
sible excited core component, the disagreement could be 
due to the spectator core approximation in the standard 
CDCC theory. Since XCDCC can handle excited core 
components, this data is re-examined. 

The breakup angular distribution data and the asso- 
ciated theory prediction for the lower energy bin (0.5- 
3.0 MeV) and the higher energy bin (3.0-5.5 MeV) are 
presented in Fi gs. EH and 2] (the equivalent of Figures 3b 
and 3c of Ref. Firstly, the CDCC calculations of 

Ref. were redone, with a higher CDCC bin density. 
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FIG. 3: n Be breakup at 60 MeV/nucleon with the relative 
energy between breakup fragments in the range 0.5-3.0 MeV. 
Experimental data are from Ref. [Til ]. 
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FIG. 4: (Color online) 11 Be breakup at 60 MeV/nucleon with 
the relative energy between breakup fragments in the range 
3.0-5.5 MeV. Experimental data are from Ref. [ill ]. 



We find that converged CDCC calculations, for the lower 
energy bin, do in fact agree well with the data (dashed 
line in Fig. [3]) , contrary to what is presented in Ref. [ll[ . 
Results do not change significantly when core excitation 
is included with XCDCC (solid line in Fig. [3]), as could be 
expected expected due to the resonance: structure in this 
energy region. One can conclude that a single-particle 
description for the first resonance is adequate. 

The data for the higher energy bin is not well described 
within the single particle CDCC model (dashed line in 
Fig. 01 . To see if this discrepancy can be explained by ex- 
cited core contributions, we include the excited 10 Be(2 + ) 
components in the reaction mechanism, within XCDCC 
(solid line). As shown in Fig.[5J core excitation lowers the 
cross section but does not significantly change the shape 
of the distribution. It becomes clear that core excitation 
does not help to reproduce the shape of the higher en- 
ergy angular distribution. The main reason for this is 
that for the 11 Be coupled channel model of [2(|, most of 
the breakup ends up in the 10 Be ground state. Fig. [4] 
also shows the breakup cross section to the + and 2 + 
states of 10 Be (red/dotted line and the dot-dashed/blue 
line respectively). We see that whereas the ground state 



distribution has the original shape of the CDCC calcu- 
lation, the shape of the distribution to the excited state 
reproduces the data (to illustrate this fact, we show the 
breakup cross section to 10 Bc(2 + ) multiplied by 10 by 
the dashed/blue line). The reason for this maybe that 
the large number of resonances in this region are not re- 
produced well by our particle-rotor model for 11 Be. The 
only resonance that appears in this model is the 3/2 + , for 
which the width is not narrow enough to attract signifi- 
cant cross section. Some suggest that more exotic struc- 
tures are responsible for resonances in this region [24| . 
Without exotic resonances built on excited core compo- 
nents in our structure model, the breakup cross section 
is still dominated by ground state 10 Be fragments. 



IV. CONCLUDING REMARKS 

A consistent analysis of reactions involving the halo 
nucleus n Be on protons, at two intermediate energies 
(~40 and ~60 MeV/nucleon) are performed and com- 
pared with data. An optical model approach, based on 
a cluster folding potential constructed from the 10 Be+p 
potential fitted to the appropriate elastic data, is un- 
able to describe the 11 Be elastic data. The inclusion of 
breakup effects improve the description, but theoretical 
predictions still overestimate the elastic cross section at 
larger scattering angles. The inclusion of core excitation 
does not affect the elastic distribution significantly. Note 
however that these results include no artificial renormal- 
ization of the optical potential. Elastic scattering ex- 
periments with radioactive beams at large facilities have 
repeatedly been undermined. The fact that the best re- 
action models are still unable to fully describe the mecha- 
nisms for the n Be case, shows the need for a more varied 
and better elastic scattering experimental program. 

In this work we also study the breakup channel explic- 
itly. Core excitation in the description of the continuum, 
within XCDCC, produces a slight modification of the dis- 
tribution. These breakup calculations are compared to 
the data at 63.7 MeV/nucleon, for two energy bins 0.5- 
3.0 MeV and 3.0-5.5 MeV. For the lower energy bin, the 
shape of the angular distribution is well reproduced by 
the models. The same cannot be said for the higher en- 
ergy bin. 

The XCDCC calculations predict breakup states to 
specific states of the core 10 Be. This level of detail is 
still not available in the data, but it could be helpful in- 
formation, even at an integrated level, to identify possible 
causes for the remaining disagreement with the data. 

Another important point is related to the basis used 
to describe the breakup states. As discussed in Ref. [25| , 
within CDCC, one can describe the three body final state 
continuum 10 Be+n+p in the 11 Be continuum basis or in 
the deuteron continuum basis. In this work we used the 
11 Be basis. Work in Ref. [25| shows that in practice the 
two choices do not provide the same result. Efforts are 
underway to tackle this problem within a Faddeev frame- 



5 



work [2^ | . These results may have important implications 
to the theory-experiment mismatch. 

Acknowledgements 

We thank the high performance computing center 
(HPCC) at MSU for the use of their facilities. This 



work is supported by NSCL, Michigan State University, 
the National Science Foundation through grant PHY- 
0555893, by NNSA through US DOE Cooperative Agree- 
ment DEFC03-03NA00143 at Rutgers University, and by 
the US DOE under contract No. DE-FG02-96ER40983 
at the University of Tennessee. 



[1] J. S. Al-Khalili and F. M. Nunes, J. Phys. G 29, R89 

(2003). [14 

[2] Y. Sakuragi, M. Yahiro, and M. Kamimura, Prog. Theor. 

Phys. Suppl. 89, 136 (1986). [15 

[3] F. M. Nunes and I. J. Thompson, Phys. Rev. C 59, 2652 

(1999). [16 

[4] J. A. Tostevin, D. Bazin, B. A. Brown, F. Glasmacher, [17 
P. G. Hansen, V. Maddalena, A. Navin, and B. M. Sher- 
rill, Phys. Rev. C 66, 024607 (2002). [18 

[5] A. M. Moro, R. Crespo, F. M. Nunes, and I. J. Thomp- 
son, Phys. Rev. C 67, 047602 (2003). [19 

[6] F. Matsumoto, E. Hiyama, K. Ogata, Y. Iseri, M. 

Kamimura, S. Chiba, and M. Yahiro, Phys. Rev. C 70, [20 
061601(R) (2004). 

[7] K. Ogata, S. Hashimoto, Y. Iseri, M. Kamimura, and M. [21 
Yahiro, Phys. Rev. C 73, 024605 (2006). [22 

[8] N. C. Summers, F. M. Nunes, and I. J. Fhompson, Phys. [23 
Rev. C 73, 031603(R) (2006). 

[9] N. C. Summers, F. M. Nunes, and I. J. Fhompson, Phys. [24 
Rev. C 74, 014606 (2006). [25 
[10] M. D. Cortina-Gil et al, Phys. Lett. B 401, 9 (1997). 
[11] A. Shrivastava et al, Phys. Lett. B 596, 54 (2004). [26] 
[12] J. S. Winfield et al., Nucl. Phys. A683, 48 (2001). 
[13] J.-P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. 



D. Childs, and Z. Vrcelj, Phys. Rev. 

and J. A. Fostevin, Phys. 



C 15, 10 (1977). 
W. W. Daehnick, J. 
C 21, 2253 (1980). 
R. C. Johnson, J. S. Al-Khalili, 
Rev. Lett. 79, 2771 (1997). 
F. Delaunay et al., in preparation. 

N. Keeley, N. Alamanos and V. Lapoux, Phys. Rev. C 
69, 064604 (2004). 

N. C. Summers, S. D. Pain et al, submitted to Phys. 
Lett. B. 

F. M. Nunes, J. A. Christley, I. J. Thompson, R. C. John- 
son and V. D. Efros, Nucl. Phys. A609, 43 (1996). 
F. M. Nunes, I. J. Fhompson, and R. C. Johnson, Nucl. 
Phys. A596, 171 (1996). 

H. Iwasaki et al, Phys. Lett. B 481, 7 (2000). 

V. Lapoux et al., submitted to Phys. Lett. B. 

R. Crespo, J. C. Fernandes, V. Lapoux, A. M. Moro, 

F. M. Nunes, private communication (2003). 

F. Cappuzzello et al, Phys. Lett. B 516, 21 (2001). 

A. M. Moro and F. M. Nunes, Nucl. Phys. A767, 138 

(2006). 

A.Deltuva, A.M. Moro, F.M. Nunes and A.C. Fonseca, 
DNP meeting, Nashville 24-27 October 2006. 



